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Abstract 

It is shown that the main variable Z of the Null Surface For- 
mulation of GR is the generating function of a constrained Lagrange 
submanifold that lives on the energy surface H = and that its level 
surfaces Z = const, are Legendre submanifolds on that energy surface. 
By globally extending the Lagrange submanifold over T*M one obtains 
a generalized generating function Z (a generating family). Thus, the 
singularity structure of the wavefronts can be obtained by studying 
the projection map of the Legendre submanifolds to the configuration 
space. The behaviour of the variable Z at the caustic points is analysed. 
It is shown that except for Minkowski space a single function Z(x, £, £) 
cannot generate the conformal structure of a radiative space-time. 

I. INTRODUCTION 

The study of null geodesies on a space-time plays an important role in general 
relativity. Gravitational lensing, either cosmological, weak or micro lensing, is the 
study of the behaviour of null rays on a curved space-time. Null geodesies are used 
to define the notion of singularities on a space-time, to define the null boundary of 
a spacetime representing compact objects, etc. 

In the last several years a formalism has been developed where null surfaces play 
a dynamical role replacing the metric as a basic variable . The goal of the Null 
Surface Formulation of GR, or NSF for short, is to introduce a new variable such 
that from its knowledge one can obtain all the conformal structure of the space-time. 
Field equations equivalent to the Einstein's equation then determine the dynamical 
evolution of those surfaces. By casting GR as a theory of surfaces rather than a 



theory of fields the NSF provides a completely new point of view with emphasis on 
the geometrical character of the theory. The basic variable is a function Z(x a , £, Q 
with x a representing points on the spacetime and ((, £) parametrizing the sphere of 
null directions. At each point on the space-time the function Z satisfies 

g ab (x)d a Z(x,(,Od b Z(x,(,Q = 0, (1) 

and the level surfaces of this function, namely Z = const, are null hypersurfaces 
on the space-time. The reader should be aware that the above construction is done 
at a local level and that in general it might not be possible to find a single func- 
tion satisfying these conditions on the whole space-time. Weyl curvature induces 
self-intersections and caustics on null congruences. Thus, even if one locally ob- 
tains a smooth hypersurface, extending such surface along the generators of its null 
geodesies will fail to be smooth. This generalized null surface is called a wavefront 
and cannot be described as the level surface of a single function Z. 

The dynamics of the NSF is imposed as field equations for Z and another scalar 
Q (the conformal factor) which are equivalent to the vacuum Einstein equations. 
The global regular solutions (in a suitable way to be specified in the conclusions) 
of those equations yield a radiative space-time, i.e., a space-time representing self 
interacting gravitational radiation. We see once again that even from a dynamical 
point of view it is very important to study the global behaviour of this variable, 
namely the behaviour of wavefronts in General Relativity. 

The purpose of this work is to provide a framework (at a kinematical level) to 
discuss global behaviour of the basic variable Z. We want to analyze the singularities 
of our variable. We want to see under what circumstances a single function Z suffices 
to construct the entire conformal structure, or how many different functions must 
be given to cover the space-time. We only consider a specific class of space-times, 
asymptotically flat space-times along future null directions. This class of space-times 
represents isolated sources that may emit gravitational radiation. 

The idea is to study the global solutions of eq.(|l|). Note that Z can be thought as 
the action of the Hamilton- Jacobi equation for the Hamiltonian H(x,p) = g ab p a Pb- 
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Since equation ([D can be written as 

H(x,d a Z) = 0, 

Z is the action of the time-independent Hamilton- Jacobi equation and the study of 
the unicity of the solution and its global properties can be carried out using the 
tools of analytical mechanics. It is worth mentioning that the study of the solutions 
of the Hamilton Jacobi equations led to the development of the theory of Lagrange 
submanifolds on cotangent bundles and the loss of unicity on the solutions is directly 
related with the singularities of the projection map of these submanifolds onto the 
configuration space 

In this paper we reintroduce our variable Z as the generating family of a con- 
strained Lagrange submanifold that lives on the energy surface H = and show that 
its level surfaces are Legendre submanifolds. Thus, the singularity structure of the 
wavefronts can be obtained by studying the projection maps to the configuration 
space. We thus, define the caustic set as the points on the Lagrange or Legen- 
dre submanifold with singular projection and the projection of those points as the 
caustics. Since Lagrange and Legendre submanifolds are smooth surfaces in T*M 
this work suggests that one can redefine our variable in a way that is free from the 
singularities and self-intersections that are naturally associated with characteristic 
wavefronts in GR. 

In Section [D] we introduce the necessary mathematical background needed for 
this work. In this context we also prove that the hypersurfaces of a constrained 
Lagrange submanifold defined as the restriction of this Lagrange submanifold to 
the level surfaces of its generating family are Legendre submanifolds on the energy 
surface H = const. 

In Section |HI| we give a specific example to clarify certain results of the previous 
section which are technically involved. We show how to construct the space-time 
wavefronts if the generating family is given and conversely, how to reconstruct the 
generating family if several smooth pieces of the wavefront are available. Note that 
we have left aside the study of the caustic points on the wavefronts and the caustic 
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set of the Legendre submanifold that generate those points. This is done in the next 
section. 



In Section [TV] we study the singularity structure of our variable Z and the main 
results are found. We show that the caustic points are obtained by choosing the 
points were 35Z, the parameter space laplacian of Z, blows up. We also show that at 
those points (Z, <5Z, <3Z) remain finite whereas d 2 Z either vanishes or diverges. Us- 
ing available singularity theorems we find as a proposition that except for Minkowski 
space a single function Z(x, £, Q cannot generate the conformal structure of a radia- 
tive space-time. Thus, in order to properly study the global behaviour of the main 
variable in the NSF one must abandon the idea of using a single function on the 
space-time and instead one has to think of our variable as a generating family Z of 
a Lagrange submanifold on the cotangent bundle of the space-time. We close this 
work with some comments of how to deal with the dynamics of the new variable. 

II. LAGRANGE AND LEGENDRE SUBMANIFOLDS. 

In this section we review the notions of Lagrange and Legendre manifolds in a 
given cotangent bundle T*M of an n-dimensional manifold M. This subject has been 
fully developed and applied to different fields in the past twenty years. Exhaustive 
treatises at a high mathematical level and/or with beautiful applications to different 
areas in physics can be found in the literature The review here presented is 

tailored to our particular needs and by no means can be considered as a substitute 
for the standard references in the field. 

In subsections A and B, we present some definitions and introduce the concept 
of constrained Lagrange submanifold in order to reinterpret, in the next section, our 
variable Z as the generating family of a Lagrange manifold. Moreover, we prove 
in proposition |11.7| that the hypersurfaces of a constrained Lagrange submanifold 



defined as the restriction of the Lagrange submanifold to the level surfaces of its 
generating family are Legendre submanifolds on the energy surface H = const. 
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A. Lagrange manifolds 



Recall that (P, uj) is a symplectic manifold if P is an even-dimensional differ- 
entiate manifold and uj is a closed nondegenerate differential 2-form on P. We 
consider a particular kind of submanifolds of P called Lagrange manifolds. 

Definition II. 1 Let (P,uj) be a symplectic manifold of dim P = 2n, a manifold L 
smoothly embedded by a map e : L — > P is called a Lagrange submanifold of P if the 
pull-back to L of the symplectic form uj on P by e vanishes on L 

e*uj = 0, 

and L is of maximal possible dimension compatible with the symplectic structure uj, 
i.e. dim L = n. 

For these submanifolds, we introduce functions called generating functions as 
follows. 

Definition II. 2 Let (P,uj) be a symplectic manifold, L a Lagrange submanifold, 
and e : L — > P an embedding. Since locally, uj = —dn, then e*uj = —d(e*K) = 0, 
so e*K = dS for a function S : L — > R (locally defined). We call S a generating 
function for L. 

From now on we will restrict ourselves to a particular class of symplectic man- 
ifolds, the cotangent bundle of an n-dimensional manifold M, denoted by T*M. 
This bundle can be assigned local coordinates (q\Pi) with (q l ) representing points 
of M and pi the local coordinates of the covectors at the point (q l ). In these local 
coordinates the closed nondegenerate differential 2-form u on T*M can be written 
as uj = dq l A dp,. 

If L is a Lagrange submanifold of (T*M, oj), then the projection map 7r : T*M — > 
M given by 7r(<f,pj) = q % , induces a map W = tx o e called the Lagrange map. The 
set of points where the rank of 7f* drops are called the singular set and the image of 
this set is called the caustic. 
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Notice also that if S : M — > R, then the graph of dS, given in local coordinates 

by 

tf,p t )eT*M:p i = -£\. (2) 



is a Lagrange submanifold, as can be easily verified. Then, S is the generating 
function of the Lagrange manifold and 7f is a diffeomorphism. 

The converse is also true, if W is locally a diffeomorphism, then L is the graph of 
dS, where S : M — > R. In this case S is only locally defined. 

Now, consider the hamiltonian system (T*M,uj,H), where H : T*M —>■ R is 
a hamiltonian function. The Lagrange submanifolds that we are interested on are 
those that can be considered as submanifolds of the energy hypersurface H defined 
by H = const. We shall refer to them as constrained Lagrange submanifolds. 

Definition II. 3 Let L be a Lagrange submanifold of T*M and H a Hamiltonian 
function, we say that L is a constrained Lagrange manifold if L C H, where H is 
an energy surface. 

Constrained Lagrange manifolds have very interesting properties. They are in- 
variant under the flow of the hamiltonian vector field Xh ( 0, Proposition 5.3.32). 
This, can be easily proved using the fact that Xh is tangent to the hypersurface 
H = const, and that L is of maximal dimension . Thus, Xh is tangent to L and the 
hamiltonian flow preserves L. 

If L is the graph of dS, where S : M — > R, then its generating function S must 
satisfy the time-independent Hamilton- Jacobi equation 

H ( gi ' w) = const ' (3) 

i.e., the generating function S = S(q l ) is the action of the Hamilton- Jacobi equation. 
Conversely, a solution of (|3|) locally defines a constrained Lagrange manifold with a 
diffeomorphic projection to M. 

However, in general (and in the problem we want to address) L will not be 
globally diffeomorphic to M. How do we handle this situation? Since L is a smooth 
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hypersurface on T*M contained in the hamiltonian flow, at most points the Lagrange 
projection will be a diffeomorphism (since the rank of W* cannot drop more than 
n — 2 on a set of points of zero measure with respect to the topology of M). We 
thus identify three different regions in L: 

1. smooth open regions of L diffeomorphic to M, and thus local graphs of a single 
function S, 

2. points of L where the normal to L is "horizontal" in T*M, i.e., the singular 
set, which divides 1 from 

3. smooth open regions of L that project down to the same open neighborhood 
of M (these points of M have more than one preimage). These regions are 
generated by several smooth functions Si. 

From the point of view of Hamilton Jacobi theory unique solutions to @ yield re- 
gions 1, multivalued solutions to @ yield regions 3 and singularities of the solutions 
yield the singular set 2. 

Let us briefly analize how to construct the three different regions from the solu- 
tions to the time independent Hamilton- Jacobi equation. 

Since the Hamiltonian flow is tangent to L, we need to solve Hamilton's canonical 
equations in order to generate the Lagrange manifold. A beautiful method to solve 
those equations is to find a generating function of a canonical transformation such 
that in the new variables the hamiltonian is independent of the new variables. 

We recall that a canonical transformation is a diffeomorphism in T*M that 
preserves the simplectic structure. Denoting the new variables by (Q l ,Pj). It is 
then easy to show that 

Pi dg i -P i dQ i = dS( 9 i ,Q»'). (4) 

i.e. the difference between the canonical 1-forms associated with the two coordinates 
is exact. The function S is called the generating function of the canonical transfor- 
mation. Note that in the context of Lagrange manifolds this function S is not the 
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generating function of a Lagrange manifold, it is called generating family, and that 
in Catastrophe Theory its normal forms are called universal unfoldings ||. Hence 
in what follows we shall refer to S as the generating family. 

If the transformation is such that the Hamiltonian becomes constant in the new 
variables, then the new Hamilton's equations are trivially solved, i.e. Pj = j3j, Q l = 
a\ where {a\(3j) are arbitrary constants. Using (H) we get 

dQ* Pl 
dS_ _ 
dq i ^ % 

and this yields the Hamiltonian flow in the variables {q l ,Pj)- 

The question is how to find this very specific generating family S. The answer 
comes from the Hamilton- Jacobi theorem. Given the differential equation 

"{*•§) =caut (5) 

the complete integral of this equation, S = S(q % ,ati), with a i: i = l...n, arbitrary 

constants, is a generating function of a canonical transformation. To see this we set 

d 2 S 



dS 

Q { = o? and P i = —. Then, if 

da; 



^ 0, Hamilton's canonical equations can 



dctidq- 

be solved by quadratures. (Jacobi's Theorem PJT0|] 
It is then clear that the set L defined by 

L= {{Q\Pj)\ Pj = Pj} 

is a Lagrange manifold for each (3j. In particular, we can choose Pj = 0, or equiva- 
lent ly 

dS(q i ,a j )_ Q 



dcti 

Then, given S(q l , a{) a complete integral of @, with ai being arbitrary constants 
for I = l...n, we can obtain L locally as the graph of dS by solving first, if it is posible, 
a>i = ai(q l ) from the equations 



and then defining S(q l ) = S(q\ ctj(q 1 )). This can be guaranteed if the rank of the 
system (|[) is r = n in the variables The Lagrange submanifold L is described 
by setting 

Pi 



ds_ 

dq i ' 



and the Lagrange map n is a diffeomorphism (i.e. the rank of W* = n). 

If r = k < n then there exists aj = ctj(q l ), for J = l...k and S = S(q\ai), 
for / = k + l...n. In this case the rank of W* drops at some points and this can be 
related to the presence of these parameters ay, for / = k + l...n. 

Then, the solution of (0), S(q\ai), defines a Lagrange submanifold L C H 
embedded into T*M by setting: 

_ dS_ 

^ % dq { 

and imposing the constraints 

dS 



1 < i < n, 



(7) 







dai 



(8) 



Since the rank of (§) is n — k in the q 1 variables, then there exist q 1 = q I {a^q J ), 
and L and 7f(L) are parametrized by (ai, q J ). The derivative 7f* can be written as 

/ dq 1 dq 1 \ 



da! dqj 
I 



where / is the identity matrix k x k, therefore it is clear that the rank of n* > k 
and it shall be strictly less than n when 

dq 1 



0. 



The set of singular points and thus the caustic set shall be isolated points, curves or 
in general a set of points of zero measure with respect to the topology of M, since 
the rank of W* cannot drop more than n — 2. This assertion can be easily understood 
since there are two vector in T*L that can not vanish under the projection, one is Xh 
and the other is the dual of dS. Notice that we can nevertheless write aj = cti(q l ) 
if we allow aj to be multivalued functions. As a consequence we obtain multivalued 
generating functions Si. 
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B. Legendre manifolds 



Odd-dimensional manifolds do not admit a symplectic structure. The analogue 
of a symplectic structure for odd- dimensional manifolds is a contact structure. 

Definition II. 4 A contact manifold is a pair (P,u)) ; consisting of an odd- 
dimensional manifold P and a closed 2-form Cj of maximal rank on this manifold. 
An exact contact manifold (P, k) consists of a (2n — 1)- dimensional manifold P and 
a 1-form k on P such that Cj = —dk is of maximal rank on P. 

Moreover, we can define a submanifold analogous to a Lagrange manifold, N of 
P called a Legendre submanifold. 

Definition II. 5 Let (P, k) be a contact manifold of dimension 2n — 1, a (n — 1)- 

dimensional manifold N such that 

e*k = 0, 

with e : iV — > P an embedding, is called a Legendre submanifold of P. 

Now, consider the hamiltonian system (T*M,u, H), the next proposition ensures 
that we can find, in a natural way, a contact submanifold of T*M. 

Proposition II. 6 Let (T*M, u, H) be a Hamiltonian system and H a regular energy 
surface, defined by H = const. Then (H,i*u) is a contact manifold, where i : H — > 
T*M is an inclusion ( Proposition 5.1.7). 

Thus the Legendre submanifolds we will consider are those that are submani- 
folds of the contact manifold (H,i*u). Moreover they are hypersurfaces of con- 
strained Lagrange manifold of a given Hamiltonian system and the projection map 
7r : T*M — > M will induce a map tt defined as tt = n o e and called Legendre map. 
The next proposition gives a description of this kind of manifolds. 

Proposition II. 7 Let L be a constrained Lagrange submanifold of the Hamiltonian 
system (T*M,u, H) and S its generating family, i.e. a solution of the time inde- 
pendent Hamilton- Jacobi equation . Then the hypersurface N of L, defined as the 
restriction of L to S = const, is a Legendre submanifold of H. 
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Proof: Given a Hamiltonian system, the Proposition IL6 ensures that (H,i*n) is 
a contact manifold and since L is a constrained Lagrange manifold, the generating 
function S(q l ,ai) of L satisfies H(q\djS) = const. Then S defines a Legendre 
submanifold N of H by setting 



dS 

Pl = — 1 < i < n, (9) 



imposing the constraints 



— dS 

S = const, and — — = 0, (10) 

OCtj 

and requiring that the rank of ([3~0|) shall be n — k + 1 in the q 1 variables. 

Recall that (|T0"| ) is an algebraic non-linear system of equations, then if we define 
the function G = (S, J^-), a solution of (0) satisfies G = 0, that is, it belongs to 
the kernel of the map G. Therefore demanding that the rank of the derivative of 
G to be n — k + 1 in the variables q 1 and in one of the q J variables, the implicit 
function Theorem guarantees that q % = q*(q J , a>i) for i E I + 1 and j £ J — 1. 

Observe that the Legendre manifold constructed in this way becomes an hyper- 
surface of L and that both are submanifold of the energy surface H.O 

The image of the Legendre map is called the wavefront and the image of the 
constrained Lagrange manifold can be considered as a wavefront family. As in the 
case of the Lagrange manifolds, the set of points where the rank of 7r* drops are 
called the singular set and the image of this set is called the caustic. If the singular 
set of L is known then intersection of this set with S = const, yields the singular 
set of the associated Legendre submanifold. 



III. EXAMPLES: A TOY MODEL 



To clarify the concepts of constrained Lagrange and Legendre submanifolds as its 
level surfaces we present in this section a specific example. Using a particular gen- 
erating family we describe the singularity structure of the wavefront. We also show 
how many pieces of smooth wavefronts are available on the configuration space for a 
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given a generating family. (Conversely, the Legendre manifold can be reconstructed 
if several pieces of smooth wavefronts are given.) We leave aside the issue of caustic 
points and caustics on the wavefronts which are analyzed in the next section. 

Since we are interested in null surfaces on a Lorentzian manifold (M, <7y) as 
Legendre submanifolds of an energy surface in T*M, we take H = -g^ (q k )piPj as 
our Hamiltonian function and consider the hypersurface H = 0. 

For simplicity, assume that n = 3 and gij = diag(l, —1,-1). Then the Hamilton- 
Jacobi equation associated with this H is 

( asy fas y _, asy = 



\dqij \dq 2 J \dq 3/ / 
i.e. the eikonal equation. The complete integral of this equation is 



where a« satisfy 



3 

S = a + } j ajqj, 
i=i 



2 2 2 n 



and a = S\ q i =0 . Hence the complete solution can be written as 



2 

« 2 



S(q l ,a!) = 5(0, a/) + aitfi + a 1 \jl - ^— J g 2 + «2<?3, 

with 7 = 1,2. The appearence of caustics shall depend on the choise of the param- 
eters a>i and 5(0, ai) or equivalently on the initial value of the integral curves of 
Xh- Choosing 5(0, ct/) = F(ati), were F are the germs of the normal forms of the 
generating functions of Lagrange manifold , we obtain the generating family for 

the singularities type A2,A 3 o A4. For example a cusp can be obtained if we choose 

4 2 

^ Oil Oil 

5(0, ai) = and «2 = • Note that the universal unfolding of a cusp is ( see 

i) 

S(a, x, y) = ±ct 4 + xa 2 + ya 

Then, the functions 
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S(qi,gs,g3,a) = + q x a + -g 2 a(4 - a 2 ) 2 + ^q 3 a 2 

92, 93, a) = -a 3 + <?i + q 3 a + ^g 2 (4 - a 2 ) ^ - ^q 2 a 2 (A - a 2 )" 2 " 

define a constrained Lagrange manifold. Note that the Taylor expansion of 
S(q\,q 2 ,q 3 ,a) around a = is precisely the universal unfolding given above. The 
explicit construction follows. 
From the equation 

|- = -a 3 + q t + q 3 a + \q 2 (A - a 2 )^ - \q 2 a 2 (A - a 2 )" 2 " = 0, (12) 
da 2 2 



we trivially obtain 



1 o, 1 1 



q 1 = a 3 - q 3 a - -g 2 (4 - a 2 ) 2 + -q 2 a 2 (A - a 2 ~) 2 - 



1 i a 2 

and by (0) we write pi = a, p 2 = -a(4 — a 2 ) 2 and p 3 = — . 

_ _ 

The map e : R 3 -> T*M 

e(q 2 ,q 3 ,a) = (a 3 - q 3 a - ^g 2 (4 - a 2 )^ + ^<? 2 « 2 (4 - a 2 )~i 

1 ^ 2n± a ^ 
a, — ck(4 — a ) 2 , — 

' 2 v ' 2 / 

is an embedding and since S satisfies ( |TTD this surface is in H = 0. The Lagrange 
map becomes 

7r'(<?2, 93, «) = (a 3 - 93« - ^92(4 - a 2 ) 2 " + ^q 2 a 2 (A - a 2 )" 2 ", g 2 , 93), (13) 
clearly in a neighbourhood of a = 

3 

^(fefea) = (a 3 - q 3 a - q 2 (l - -a 2 ),q 2 ,q 3 ), 

o 

and this yields a cusp. 

The map ([13]) is not a diffeomorphism, since the equation 

J = det(7r*) = -3a 2 + q 3 + q 2 a{4 - a 2 )~^{a 2 - 6) = 0, 

gives us the points where the map losses its rank. Solving this equation we obtain 
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ft = /(ft, a) 

= 3a 2 - q 2 a{A - a 2 )~^{a 2 - 6). 



Then the caustic set is the image of e(qz, f(q2, a), a) and the caustic is the image of 
7f J (g2, /(ft, a), a). The last is shown in FIG [p. 

Remark III.l Observe that the projection of the surface J = into configuration 
space (i.e. the caustic) divides regions on which the Lagrange map is a diffeomor- 
phism, for example the region given by a > 0, q 2 > and q% < ; from other regions 
where the Lagrange map in not injective (more than one preimage). In the first re- 
gion we may locally write a = a(q l ) from equation and the generating function 
S : M — > R. On the other hand, in the regions of non-injectivity, we obtain more 
than one function a, as solution of the equation ftUfy which in turn implies that we 
get several functions Si : M — > R. 



FIG. 1. The caustic 

The corresponding Legendre submanifold is then constructed as a level surface 

of the Lagrange manifold presented above. 

1 i a 2 

As before we write pi = a, p 2 = -a(4 — a 2 ) 2 and p% = — . From the system 
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„ Ot 1 2 1 , 1 

S = — — + q x oi + -qzat + -q 2 a{A - a ) 2 = 

^- = -a 3 + q 1 + q 3 a + \q 2 (A - a 2 )^ - \q 2 a 2 (A - a 2 )~^ = 0, 
oa 2 2 

we obtain 

1 1 3 1 

qi = -a 3 - 2g 2 (4 - a 2 )~^ , q 3 = -a 2 + q 2 a{4 - a 2 )'^ 

The map e : R 2 -> T*M 

e(q 2 ,a) = (^a 3 - 2g 2 (4 - a 2 ) _ 5,g 2 , ^a 2 + <?2«(4 - ct 2 ) - ^ 

1 2n± ^ 

a,-a(4-a ) 2 ,yJ 

is an embedding and it defines a Legendre submanifold. The Legendre map 

1 1 3 1 

7t%q 2 , a) = {-a 3 - 2g 2 (4 - a 2 )-^,q 2 , -a 2 + g 2 a(4 - a 2 )^) 

describes a wavefront, i.e. n(N). It is shown in FIG |2|a. 




FIG. 2. (a) The wavefront and (b) J = 



The derivative of the Legendre map is calculated as 

( -af(a,q 2 ) -2{A-a 2 )-\\ 
1 



7T 



\ 2/(a,g 2 ) a(4-a 2 ) 3 ^ 
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4g 2 + 3a(4 — a 2 )a 

where f(a, q 2 ) = 3 . Clearly, the map % is not a diffeomorphism 

2(4 — s 2 ) 2 

since the rank of its differential drops when 

4g 2 + 3a(4-a 2 )5 =0, 

i.e. on the curve (a,q 2 (a)), a; G [—2,2]. This region (the caustic set) is depicted 
in FIG 0b. The caustic defined by these points (i.e. the image of the curve given 
above under tt) is given by 

3ct 3 3ct 2 

qi = - a ( a 2 -3), q 2 = — ^(4-a 2 )t and q 3 = -^-(a 2 -2), ae[-2,2]. 

and it is drawn in FIG |3|a. 



FIG. 3. (a) The caustic and (b) the null geodesies 

This caustic, together with some integral curves of the non-vanishing null vector 

l l = = f-2(4-a 2 )"^,l,a(4-a 2 )"^ are shown in FIG Bb. Note that the 

caustic is the envelope of the null geodesies with tangent vector l l and that the 

tangent vectors to the curves given q 2 = const vanish on the caustic. The vanishing 

of these tangent vectors M % = — — can be seen in FIG HI, where the caustic and 

da 

some q 2 = const curves (surfaces in higher dimensions) are shown. 
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FIG. 4. The curves q2 = const 



Observe that in the region given by —2 < a < 2 and g 2 > 4 the Legendre map is 
a diffeomorphism (see FIG 0). Thus, for a sufficiently large g 2 , the function S\ q2=const 
defines a submanifold Q diffeomorphic to Q ~ R x [—A, A], being A a positive and 
large constant (see FIG ||) . 




FIG. 5. A diffeomorphic region 
From this example we conclude that: 
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Remark III. 2 Since there are regions outside the caustic set where the generating 
function S(qi) is multiple-valued, we have several pieces of a smooth wavefront. Con- 
versely if we have several smooth pieces of a wavefront, using ^), we can reconstruct 
the Legendre manifold except at the caustic points. 



IV. THE FUNCTION Z IN ASYMPTOTICALLY FLAT SPACE- TIME 

For simplicity we will consider asymptotically flat space times with a null bound- 



ary that represents the end points of future directed null geodesies [|TTJ. Those space 
times represent compact objects that can emit gravitational radiation. 

To define our variable Z we consider the intersection of the future null cone from 
x a with the null boundary X + . Introducing Bondi coordinates (u, Q on 1 + (with 
u representing a Killing time and (£, Q being stereographic coordinates on the unit 
sphere) this intersection is locally given by the equation 

u = Z(x a ,(X). (14) 

Thus, for fixed values of x a , the function Z yields the parametric description of the 
light cone cuts of I + . 

On Minkowski spacetime the I.e. cuts adopt a very simple form, u = x a l a with 
l a a null covector constructed from the spherical harmonics Y , Yx-i, Y 10 , Y\j_. 

In general a cut is a complicated surface with caustics, self-intersections, etc. 
However, it can be shown that for regular space-times the index number is always 
one. Therefore, a I.e. cut (as complicated as it might be) is always a continuous 
deformation of the sphere of null directions above each point x a . It can also be 
shown that generically the I.e. cuts can only have two kinds of singularities, cusps 
and swallowtails since they represent the projection of 2-dim Legendre submanifolds 
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A second meaning can be assigned to our variable Z(x a ,(,0- Fixing a point 
(u, C, C) °f the collection of interior points x a that satisfy 

Z(x a , C, C) = u = const, (15) 



form the past null cone of (u, C, C)- Moreover, from knowledge of Z we can construct 
a null coordinate system as follows. 

Starting with our variable and taking (£, () derivatives of Z(x a , (, Q we construct 
the following set of scalars, 



For fixed values of (£, C) they define a coordinate system with the following geometric 
meaning, 

• u = const, denotes the past null cone from (u, £, £). 

• (u>,u) = const, single out a null geodesic on that surface. 

• R = const, identifies a point on that geodesic. 

However, one knows that null cones can develop caustics and singularities. One 
also knows that past those singularities the null cone is no longer smooth (it is called 
a wavefront) and thus, a null coordinate system like the one above breaks down past 
those singular points. Since the main goal of the NSF is to replace the metric with 
a function Z such that its level surfaces are past null cones from X + , we imme- 
diately face a non trivial problem: if the null cones develop self-intersections and 
singularities that cannot be analyzed with a single function Z then the construction 
given above is only valid on a neighborhood of X + . However, we also know that null 
wave fronts are projections of Legendre submanifolds that live on T*(M). It would 
then appear that the best way to deal with this lack of smoothness is to think of 
our variable as the generating family Z of a constrained lagrange submanifold. An 
outline of this construction is presented below. 

As was done before we assume that Z is a solution of the equation 



9 i (x a ,CX) = (O°,9 + ,0-,9 1 ) = (u,u,Q,R) = (Z,QZ,QZ,QQZ). 



(16) 



H(x a ,d b Z)=g ab Z, a Z, b =0. 



(17) 



with 



g ab a metric that is asymptotically flat. 
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In a neighborhood of X + the solution to this equation yields a single function 
Z and its level surfaces Z = const, describe the past null cones from points at the 
null boundary. This follows from the fact that the (unphysical) metric near X + is 
"almost" conformally flat and thus the past null cones are free from caustics and 
singularities. 

Since Z(x a , Co, Co) is a smooth function on this region, we can choose S = 
Z(x a , Co, Co) as the generating function of a constrained Lagrange manifold L 

N=Ux a )Vb = ^):e*K = dz\. (18) 



Note that the manifold described above is equivalent to (|2]) and the Proposition |II.7 



ensures that the surface N defined by Z = const is a Legendre submanifold of the 
energy surface given by H = 0, i.e. 

N = j(x a ,p 6 = 0) : e*k = d(e*Z) = oj . 

Thus, we have constructed a constrained Legendre submanifold N and a con- 
strained Lagrange submanifold L of H using our fundamental variable Z(x a , Co, Co)- 

The idea is to extend this construction to regions were caustics develop. As we 
mentioned before, in these regions the Lagrange submanifold is not diffeomorphic 
to its projection. We will thus assume that a solution of H(x a , dbZ) = can be 
written as Z = Z(x a , Co, Co, a i) with 1—1,2. The function Z depends, at most, on 
two parameters since the rank of the projection map cannot drop more than two [Q. 
The constrained Lagrange submanifold L as described in section [11 A| is given by 

dZ 

Pa 



dx° 



together with the constraint 



Observe that if we can solve (O) uniquely for ai, i.e. oti = ai(x b ), then Z = Z and 
we are back in the previous diffeomorphic region. In general, one will obtain multi- 
valued solutions of (|i9|). Inserting the different solutions of ai into Z(x a , Co, Co, 
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one obtains a multiple- valued function Z(x a , £ , Co)- The Legendre submanifold is 
obtained by setting Z = const. Conversely, if several functions Zi are given, one can 
reconstruct the Lagrange submanifold by imposing ( |18|) on the different Z's. The 
construction defines the lagrange submanifold except for the caustic set. 

Finally, we would like to determine under what circumstances it is possible to find 
a single function Z that would yield for us a global coordinate system (u, R,w,w) 
on an asymptotically flat space time. In other words we want to know if there exists 
space times that are diffeomorphic to the corresponding Lagrange manifolds. At the 
same time we would like to know when and how this coordinate system breaks down 
due to the presence of conjugate points. We are therefore interested in describing 
the relationship between our fundamental variable Z and the loss of the rank of the 
derivative of the Legendre map n, i.e. we want to described the singular set in terms 
of Z. 

When the Lagrange manifold is a constrained one, the loss of rank of the Lagrange 
map indicates the non-existence of global solutions of the Hamilton- Jacobi equation 
and the loss of rank of the associated Legendre map is related to the existence of 
conjugate points of a congruence of null geodesies. 

In order to clarify this assertion, we consider the local description of the wavefront 
(the projection of the Legendre manifold). We assume that the wavefront is locally 
described by 

x a = f a (u ,s,w,wXoXo)- 

with s an affine length. The vectors 

df a df a - df a 

L a = M a = -^— and M a = 

as ow ow 

are tangent to the wavefront. L a is directed along the null geodesies whereas M a 

and M a are geodesic deviation vectors. 

The derivative of the Legendre map losses its rank when these three vectors 

become linearly dependent. This dependence is related to the existence of conjugate 

points on the congruence of null geodesic with apex at X + and null tangent vector 

L a as follows. 
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We introduce the parallelly propagated null triad {l a , m a , m a } , satisfying 



l a m a 

m a m a 

l a V„m b 



0. 



(20) 



In terms of this triad 



therefore this set of vectors becomes linearly dependent when 



(21) 



(22) 



On the other hand, this quantity is related to the divergence p, and the shear a of 
the congruence with apex in X + . To see this , consider the optical parameters 



p = m a m V Jb and a = m a m VJb- 

Using eq. ([H]) together with the fact that M a is Lie propagated along the null direc- 
tion L a we get ( 0) 



a = —D - 
A \r] 



P 



DA 
~2A 



with A = (££ — rjfj) and where we have used the fact that p is real. 

Hence, at the points where the Legendre map looses its rank, the divergence 
of the congruence becomes unbounded, i.e. limp s ^ SQ = oo, where s is the affine 
parameter and Sq corresponds to a conjugate point. 

Summarizing, the loss of rank of the Legendre map, i.e. the development of a 
caustic is directly related to the existence of a conjugate point on the null congruence 
with apex at I + . 

Since Z is a unique function near X + , (u, R,w,w) is a well behaved coordinate 
system in that region. The question is what happens to our coordinates as we 
approach a generic conjugate point. 
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Lemma IV. 1 Let Z(x a ,(,C) be our basic variable near X + such that Z = de- 
scribes the past null cone from (u,(X)- Then at the first conjugate point of this null 
congruence the scalar R = 9SZ goes to — oo. 

Proof: Given a null geodesic labeled by (u, w, w, C, C) we introduce an affine length 
s and two null congruences that contain this geodesic; 1) the future null cone with 
apex at So and 2) the past null cone with apex at (u, C,C)- The divergence of the 
first congruence at 1 + is related to the value of the scalar R(sq) in the following way 
ITS]] . The divergence of the generating vectors of the first cone is defined by 

Pi = m a m h V a F b . 

where F = F(Q, u, x l ) and F = describes the null cone. Near X the function F 
can be written as 

F = F°{u, x l ) + VLF\u, x') + 0{Q 2 ). 
where F° = u — Z. Then 

Pi(s ,J + ) =m a m b V a (u- Z) 

= p B -R{s ). (23) 

where ps is the divergence of a Bondi congruence. 

On the other hand, the Sachs- Penrose reciprocity theorem, states that: 

Theorem: Assume that X% and X2 are two matrices whose elements are the tetrad 
components of the two complex deviation vectors associated with a null cone congru- 
ence with apex at a point p\ and P2 respectively, then 

Xi( atp 2 ) = -X 2 ( atpi). 

Hence using this theorem, we may assert that if the past null cone from I + has 
a conjugate point at an affine distance s = s* then the future null cone from s* has 
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a conjugate point at X + . Then in the limit when sq — > s*, we find that pi(so,X + ) 
goes to +00 and from(|2"5|) R(s*) = — 00 □ 

The first consequence of this lemma is that our coordinate system is well defined 
in the domain R G (— 00, 00), in other words from our coordinate system we can not 
detect the caustics that arise in the past null cones as we move into the space-time 

The lemma is also useful to answer the question we possed before, namely, if there 
exist asymptotically flat space times that can be covered with a global canonical 
coordinate system constructed from a single function Z. Using proposition 4.4.5 



TT| which states that any null congruence along a geodesic such that the afline 
length can be extended arbitrarily has a pair of conjugate points and lemma [IV. 1| 
we conclude that the coordinate system (u,w,w,R) derived from Z cannot cover a 
space-time except for Minkowski space. We conclude that 

Proposition IV. 2 A single function Z(x, £, Q cannot generate the conformal 
structure of a radiative space-time except for the Minkowski space. 

Is is easy to show that <SZ and BZ remain finite at a conjugate point. This 
follows from the fact that both are constants along a null geodesic. 

It is also of interest to analize the behaviour of the conformal factor Q, and 
A = <S 2 Z near a caustic since they generate the underlying metric of the space-time. 

In a similar way than the Lemma above, we consider the future null cone with 
apex at Sq. The Sachs Theorem tells us that at X + 

ai(s ,X + ) = m a m b VJ b = a B - A(s ). (24) 

where Ob is the shear of a Bondi congruence. 

If the future null cone of Sq has a conjugate point at X + , then the shear o~\ (sq,1 + ) 
is either or 00. It then follows from the Sachs-Penrose reciprocity theorem and 
eq. ([24] ) that |A(sq)| is either or 00 if the past null cone from X + has a conjugate 



point at sq. 

It remains to consider the conformal factor 

dR 



n 2 = g 01 :=g ab Z, a mZ. 
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b ds 



Since R(s) diverges as s approaches a conjugate point while the affine length is a 
smooth non- vanishing function along the null geodesic it follows that g 01 also blows 
up at that point. 

V. CONCLUSIONS 

We have shown that our main variable can be regarded as the generating family 
Z of a constrained Lagrange submanifold and that its level surfaces are constrained 
Legendre submanifolds that project down to past null cones from 1 + . Further- 
more, we have demonstrated that for a generic space-time this Lagrange subman- 
ifold starts diffeomorphic to the configuration space but it develops caustic sets, 
points on T*(M) where its projections are caustic points on M. Thus, except for 
Minkowski space, a single function Z on configuration space does not give the con- 
formal structure of the space-time. At the caustic points QdZ diverges. This means 
that the coordinate system constructed on the null cones is only locally defined but 
on the other hand one never sees the caustics since they are pushed out to R = — oo. 

Although the entire treatment so far has been kinematical we would like to think 
of our variable as coming from the solution of a set of field equations given on the 
space-time ||[| . 

It is clear from the previous results that the solution of those field equations 
must have multiple valuedness in order to generate the multiple branches needed 
to construct the generating family Z of the Lagrange submanifold. These solutions 
are defined in a 6-dimensional space, 4-spacetime coordinates and 2 parameters on 
the sphere, (CO- We demand the solution to be globally defined respect to the 
parameters ((, C), that is, it shall be a piecewise smooth function on the sphere (it 
could be multiple valued but always finite on the sphere). In this sence we say that 
the solutions shall be globally regular on the space of parameters. 

Alternatively, we could try to find field equations given on T*M. In this case 
the solution would yield a global generating family Z of a constrained Lagrange 
submanifold that coincides with Z in a neighbourhood of X + . This last approach 
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will be further explored. 
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